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| The thermodynamical behaviours of ferromagnetic systems in equilibrium are well studied, 
f-r^ However, the ferromagnetic system, far from equilibrium, became an interesting field of research 
in last few decades. The ferromagnetic systems in the presence of a steady magnetic field are 
r- 1 ! also studied by using standard tools of equilibrium statistical physics. The ferromagnet in 
O ■ the presence of time-dependent magnetic field, shows various interesting phenomena, explored 
rH ■ recently. An usual response of a ferromagnet in the presence of a sinusoidally oscillating magnetic 
i ! field is the hysteresis. Apart from this hysteretic response, the nonequilibrium dynamic phase 
^ transition is a very interesting phenomenon. In this chapter, the nonequilibrium dynamic phase 

■ transitions, in model ferromagnetic systems in presence of time-dependent magnetic field, are 
discussed. For this kind of nonequilibrium phase transition one cannot employ the standard 
techniques of equilibrium statistical mechanics. The recent developments in this direction are 

^ ■ mainly based on numerical simulation (Monte Carlo). The Monte Carlo simulation, of kinetic 
'"O ■ Ising model in presence of sinusoidally oscillating (in time but uniform over space) magnetic 
^ \ field, is extensively performed to study the nonequilibrium dynamic phase transition. The 
q ■ temperature variations of dynamic order parameter, dynamic specific heat, dynamic relaxation 
time etc. near the transition point are discussed. The appearance of a dynamic length scale and 
a dynamic time scale and their behaviours near the transition point are also discussed. All these 
studies indicate that this proposed dynamic transition is a nonequilibrium thermodynamic phase 
transition. The disorder (quenched) induced zero temperature (athermal) dynamic transition is 
i— I studied in random field Ising ferromagnet. The dynamic transition in the Heisenberg ferromagnet 
0^ | is also studied. The nature of this transition in the Heisenberg ferromagnet depends on the 
jy-j anisotropy and the polarisation of the applied time varying magnetic field. The anisotropic 
, Heisenberg ferromagnet in the presence of elliptically polarised magnetic field shows multiple 
dynamic transitions. This multiple dynamic transitions in anisotropic Heisenberg ferromagnet 
are discussed here. Recent experimental evidences of dynamic transitions are also discussed very 
briefly. 
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- - 1. Introduction: 

The ferromagnetic system in statistical equilibrium gives rise to well-known ferro-para phase 
transition. However, the ferromagnetic system, in the presence of a time varying external mag- 
netic field, remains far from thermodynamical/statistical equilibrium. This type of system 
became an interesting object of research over the last two decades. Particularly, if the time 
varying externally applied magnetic field is sinusoidally oscillating (in time but uniform over the 
space), it yields two major responses of the ferromagnetic systems, i.e. (i) dynamic hysteresis (ii) 
dynamic phase transition. The dynamic hysteretic responses of Ising-like systems are already 
discussed in a recent review article [1]. But recently, the another response, i.e., the dynamic 
transition has drawn much attention of researchers. It has several interesting and important (in 
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the context of nonequilibrium phase transition) aspects and can provide a simple example of 
nonequilibrium transition. In the earlier review [1], the dynamic transition was also discussed 
in brief. However, in last five years a considerable amount of research was done on it (in the 
Ising model) which was not reviewed yet. Very recently, the study of the dynamic transitions 
in the Heisenberg model gives several interesting features which are not possible to observe in 
the Ising models. Keeping this in mind, I have tried to review the recent developements in this 
field. 

This review is organised as follows: in the next section a short introduction is written 
about what a dynamic transition is and when it should be called true dynamic transition. In 
third section, the dynamic transition in the Ising ferromagnet is discussed. Different aspects of 
dynamic transitions (mainly its thermodynamic natures) and its relationship with the stochastic 
resonance are reviewed. In this review the analytic formulation of dynamic transition and the 
athermal dynamic transition. In the fourth section, the dynamic transition in the Heisenberg 
ferromagnet is discussed. This section gives a review of very recent (last 1-2 years) works on 
classical vector spin models. The results available so far, based on the experimental observations 
reviewed in the fifth section. The chapter ends with a summary given in the sixth section. 



2. What is dynamic transition ? 

If the temperature of a ferromagnetic sample increases (in absence of external mag- 
netic field), the sponteneous magnetisation vanishes at a particular temperature (Curie 
temperature) and the transition occurs from an asymmetric (ferromagnetic) phase to a 
symmetric (paramagnetic) phase. This is a very well known phenomenon and the equilib- 
rium symmetry breaking ferro-para transition is also well studied. But any nonequilibrium 
transition of any ferromagnetic sample in the presence of time-dependent magnetic field 
was unexplored even before 1990!! Tome and Oliviera [2] first observed a prototype of 
nonequilibrium dynamic transition in the numerical solution (by fourth order Runge- 
Kutta method) of mean field equation of motion 
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of the classical one component (m) ferromagnetic model in presence of a magnetic field 
varying sinusoidally (h(t) = h cos(ut)) in time. Here, the time-averaged magnetisation 
Q = § m(t)dt/r (r = 2-k/lu is the time period of the oscillating field) over a full cycle of the 
oscillating field plays the role of order parameter for this type of proposed nonequilibrium 
dynamic transition. Tome and Oliveira [2] found that Q becomes zero (dynamically dis- 
ordered phase) from a nonzero (dynamically ordered phase) value at a finite temperature 
T (which also depends on the value of amplitude h of oscillating field). They sketched 
a dynamic phase boundary in h — T plane. For the higher values of ho the transition 
was found to be of discontinuous type and that becomes continuous for lower values of 
h . They [2] also located a tricitical point (TCP) on the phase boundary which separates 
the nature (discontinuous/continuous) of these transitions. 

This phenomenon can be explained by considering the system initially kept in one 
well of a Landau type double well potential. Depending on the temperature, a definite 
amount of magnetic field is necessary to bring the system from one well to another. If the 
amplitude of the applied oscillating magnetic field is less than the required amount, the 
system oscillates in one well (where it was initially). In this situation, the magnetisation 
does not change its sign. As a result, the magnitude of the time-averaged magnetisation 
is nonzero (Q ^ 0). If one sees the plot of m(t) — h(t), it is asymmetric about the 
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m(t) = line. This gives rise to a dynamically ordered and asymmetric phase. As the 
temperature increases, the height of the barrier between two wells decreases and to push 
the system from one well to another, less amount of magnetic field is necessary. As a 
result, the magnetisation can change its sign for this amount of field. Consequently, the 
time-averaged magnetisation over a full cycle becomes zero (Q = 0) when m(t) — h(t) 
plot is symmetric. So, we get the symmetry broken dynamic phase transition, where 
the transition temperature decreases as the value of magnetic field increases. This can 
qualitatively explain the phase boundary of dynamic phase transition. Fig. 1 illustrates 
this symmetry breaking associated with the dynamic transition. 

But, can we really call this transition a true dynamic transition ? A true dynamic 
transition will be such that this transition should disappear in the static limit !! Let 
us examine logically, what happens in the static (infinitesimally small frequency) limit. 
One can simulate this situation by varying the magnetic field with infinitesimally small 
frequency. It should be noted here (in eqn.l) that the dynamics of the one component 
magnetisation (m) is purely deterministic. So, at a particular temperature, to bring the 
system from one well to another, a definite amount of field should be applied, irrespective 
of the rate of achieving the required amount. The system will wait untill it gets the 
required value of the field. If the frequency is very small, it will wait for very long time. 
But unless it gets the required value of field it will not go to another well. There is no 
stochasticity involved or any noise is present in the equation of motion which can push 
the system towards another well, irrespective of value of the applied field. So, obviously, 
in the zero frequency limit a so called dynamic transition will be observed in the case of 
meanfield study. Just by this argument one can immediately conclude that the kind of 
transition described above cannot be truely dynamic in nature!! 

After realising this, the researchers, interested in this field, tried to study the true 
dynamic transition in ferromagnetic model systems incorporating the thermal fluctuations 
as the source of noise or stochasticity, keeping in mind that this stochasticity will help 
to push the system into another well. Rao, Krishnamurthy and Pandit [3] and Dhar 
and Thomas [4] tried to observe this perfectly dynamic transition in iV-vector model 
in the N — > oo limit. They [3, 4] also tried to observe this in the kinetic Ising model 
by Monte Carlo (MC) simulation. At the same time, Lo and Pelcovits [5] studied the 
kinetic Ising model in the presence of sinusoidally oscillating magnetic field by Monte 
Carlo simulation. But, unfortunately both failed to observe the dynamic transition and 
to draw the phase boundary. Being motivated from these studies; extensive MC studies 
[4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15] were performed in the kinetic Ising model in oscillating 
magnetic field and the recent developments in this direction are discussed in details in 
the next section. It should be mentioned here that the dynamic response at very low 
frequency, essentially originates from traditional nucleation problem [16]. However, to see 
this effect one should use a frequency which is much smaller than those usually employed 
in this field of research [1]. 

Although the numerical solution of mean-field equation (Eqn. 1) cannot provide for a 
true dynamic transition, it has one importance. The linearised equation is exactly solvable 
and can be used to get some qualitative features of dynamic transition analytically. 

3. Dynamic phase transitions in the Ising ferromagnet: 

(a) Model and simulation 

To investigate the true dynamic phase transitions (incorporating the fluctuations) one 
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simple choice may be the kinetic Ising model in presence of sinusoidally oscillating (in 
time but uniform over the space) magnetic field studied by Monte Carlo simulation. For 
this choice, one can take the following Hamiltonian, 



where, the first term represents the spin-spin the Ising type ferromagnetic interaction. 
Si(±l) is Ising spin at i-th site of the lattice, J(> 0) is nearest neighbour ferromagnetic in- 
teraction strength. The second term represents the spin-field interaction. h(t) = h cos(ujt) 
is applied oscillating magnetic field. h is the amplitude and uj{— 2nf) is the angular fre- 
quency of the oscillating field. The boundary condition is periodic in all sides of the 
lattice. 

Due to the presence of second term in the Hamiltonian, the system always remain 
far from equilibrium. However, the dynamical evolution of the system can be studied by 
Metropolis algorithm [17]. 

The studies on Ising systems in the presence of an oscillating magnetic field have 
been made mostly by employing the Monte Carlo method using Metropolis single spin 
flip dynamics [17]. Starting from a (random or perfectly ordered) initial configuration of 
spins, the spin state Si(t) at any site i and in any time t for a fixed temperature T is 
updated (sequentially or randomly) with the following probability function [17] 



where AH = 2Si[T l jSj(t)+h(t)], is the change in energy due to spin flip, Kb is Boltzmann's 
constant and T is the temperature. First the instantaneous response magnetisation per 
lattice site at time t is easily calculated: m(t) = YliSi(t)/N, where N is the total number 
of spins in the lattice (N = L d if one considers a d- dimensional hypercubic lattice of 
linear size L). N such spin updates is defined as one MC step per spin(MCSS). This 
is the unit of time in this simulational study. After that, one can study the dynamical 
response of the system. In this chapter, only the dynamic transitions, out of various kinds 
of dynamical responses, will be discussed. The dynamical order parameter (as defined in 
ref[2]) can be calculated as Q = §m{t)dt/r, where r = (1//) (/ is linear frequency) is 
the time period of the applied oscillating magnetic field. The magnitude of the magnetic 
field is measured in the unit of J and the temperature is measured in the unit of J / K B . It 
should be mentioned here that the results are independent of the dynamics (Metropolis, 
Glauber etc.) employed to study the dynamic transitions. 

(b) Dynamic transition: 

The value of stabilised Q is calculated [6] for fixed values of T, ho and /. It was 
observed that for a fixed frequency /, the dynamic order parameter Q is nonzero for 
lower values of T and ho and it would vanish for higher values of T and ho. For lower 
values of temperatures and field amplitudes the system is dynamically ordered and it 
loses its dynamical ordering for higher values of T and h - This dynamical order- disorder 
transition is associated with the breaking of symmetries of to — /i loop . For lower values 
of T and ho, the magnetisation oscillates asymmetrically (around to = line) which gives 
rise to asymmetric m — h loop (resides asymmetrically in to — h plane) (see Fig.l). In 
this case Q ^ 0, and this phase is a symmetry-broken phase. When T and h become 
large, the magnetisation oscillates symmetrically (around to = line). This gives rise to 
symmetric to — h loop (resides symmetrically in to — h plane) and consequently Q — 0. 




(2) 




(3) 



4 



This phase is called symmetric phase. In this regard, one may call that this dynamic 
transition is associated with a dynamical symmetry breaking. 

Extensive Monte Carlo (MC) simulation was performed [1] to study this dynamic 
transition. MC simulation was done in the Ising ferromagnet in two and three dimen- 
sional hypercubic lattices. The dynamic transition was observed and the dynamic phase 
boundary was drawn in h — T plane taking frequency as a parameter. The dynamic 
phase boundary in both the dimensions (two and three) are qualitatively similar in na- 
ture. Along the phase boundary the dynamic transition occurs at higher temperatures for 
lower values of field amplitudes. The high temperature (low field) transitions are contin- 
uous and the low temperature (high field) transitions are discontinuous. A point (marked 
TCP in figure 2) separating these two types of transitions is called a tricritical point. An- 
other important thing should be noted in this dynamic phase boundary is, the variation 
of the phase boundary with respect to the frequency of the driving field. As Fig. 2 shows, 
the phase boundary shrinks inward as the frequency decreases. As an extrapolation of 
this senario, one may conclude that, in the static limit (zero frequency limit) the dynamic 
transition disappears. This is the important significance of true dynamic transition. 

The dynamic transition is, in fact, a manifestation of the coercivity property (one of 
the important features of hysteresis) [12, 16]. In the Q ^ phase, the m — h loop is not 
symmetric about the field axis and lies in the upper (or lower) half of the m — h plane 
depending on their initial magnetisation. A minimum magnitude of external and oppo- 
site magnetic field (coercive field) is required to change the sign of the magnetisation for 
complete reversal within the time period of the oscillating magnetic field. This magni- 
tude of the coercive field depends on the temperature T. The magnitude of coercive field 
increases as the temperature decreases. In the case of sinusoidally oscillating magnetic 
field, for a transition to a Q = phase from Q ^ phase, the field amplitude should be 
at least of the order of coercive field depending upon the temperature T. So, in a sense, 
the dynamic phase boundary (in the low frequency limit) is the coercive field variation 
with respect to the temperature. Since, h(t) = h cos(ut) and \h(t)\ < h , the phase 
boundary, is the upper bound of the coercive field variation with respect to the tempera- 
ture T. However, the difference of this upper bound increases with increasing frequency, 
because of the effectice relaxational lag (r e //, time lag of magnetisation with respect to 
the applied field). The effective relaxational lag of the magnetisation arises due to the 
intrinsic relaxation time of the system. In fact, the tricritical point Tj CP (h ,uj) on the 
phase boundary appears because of the system's failure to relax within the time period 
2h/uj of the oscillating field. The intrinsic relaxation time in the ferromagnetic phase 
decreases with lowering of temperature and below Tj CP (ho,uj), T e ff < 2-k/uo (equality at 
T = Tj CP ). So, that the magnetisation changes sign (from +m to —m) abruptly and 
consequently Q changes from a value very near to unity to zero discontinuously. This 
indicates the TCP should decrease with higher frequency. At zero temperature, the tran- 
sition is completely mechanical (purely field driven; without any thermal fluctuation) and 
can only be a discontinuous one. Above Tj CP (ho,u), the thermal fluctuations win-over 
and determine the continuous nature of the transition. There is a controversy regarding 
the existence of TCP. The details of studies on TCP will be discussed in the subsection 
where the relation between dynamic transition and stochastic resonance, is discussed. 

It should be mentioned here that another type of dynamic phase transition was studied 
[18] in the kinetic Ising model with negative pulsed magnetic field of finite duration. But 
in this chapter, the dynamic phase transition, only due to sinusoidally oscillating magnetic 
field, will be discussed. 
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(c) Is dynamic transition a phase transition ? 

After getting the dynamic transition and the dynamic phase boundary, the following 
questions naturally arise: Is dynamic transition a phase transition ? Is there any evidence 
of the divergence of time scale and length scale at the transition point ? What will be the 
behaviour of 'dynamic specific heat' and fluctuations of dynamic order parameter near the 
transition point ? The studies related to the above questions will be discussed in this 
subsection. 

The critical slowing down is an important phenomenon observed in equilibrium tran- 
sition (ferro-para) which indicates the divergence of time scale (relaxation time) at tran- 
sition point. In the MC simulation if the initial condition is all spins directed upward and 
the sinusoidal magnetic field is driving the system, the m — h loop gets stabilised after 
a transient behaviour. As a result, the dynamic order parameter Q has also a transient 
behaviour. 

It has been observed carefully [7] that the dynamic order parameter Q does not acquire 
the stable value within the first cycle of the oscillating field. It takes several cycles (of the 
oscillating field) to get stabilised i.e., it shows 'relaxation' behaviour. Starting from the 
initial (all spins are up) configuration, the Q has been calculated for various numbers (say 
n-th) of cycles of the oscillating magnetic field and plotted (inset of Fig. 3) against the 
number of cycles (n). Each value of Q shown here has been obtained by averaging over 
100 random Monte Carlo samples. Inset of Fig.3 shows a typical 'relaxation' behaviour of 
the dynamic order parameter Q. This has been plotted for fixed values of u = 2ir x 0.04, 
ho— 1.0 and T— 1.5. It shows that the dynamic order parameter Q is relaxing as the time 
(number of cycles) goes on. The best fit curve shows that the 'relaxation' is exponential 
type. So, one can write Q ~ Q Q exp(— n/T), where T is the 'relaxation' time which 
provides the 'time scale' for this prototype of nonequilibrium problem. The physical 
interpretation of T is, the number of cycles required, so that Q becomes 1/e times of its 
initial value (value of Q at starting cycle). From the exponential fitting, the 'relaxation' 
time (T) has been measured. The temperature (T) variation, for fixed values of u> and 
ho, of this 'relaxation' time T has been studied (in the disordered region of dynamic 
transition) and displayed in Fig.3. The temperature (T) variation of T are shown (Fig.3) 
for two different values of ho(— 0.5 and 1.0) and for a fixed value of uo = 2n x 0.04 here. 
From the figure (Fig.3) it is clear that the relaxation time T diverges near the dynamic 
transition point (where Q vanishes) in the both cases (h = 0.5 and 1.0). This study [7] 
first indicated that this dynamic transition is associated with a diverging 'time scale'. 

An analytical formulation of this critical slowing down (of Q) can be done [7] by solving 
the linearised mean-field equation of motion for the average magnetisation . In the limit 
of h — > and T > 1, the equation (eqn. 1) can be linearised (i.e., linearising tanh term) 
as 



where m and are two constants. The value of the dynamic order parameter Q at n— th 
cycle of the oscillating field is 




where e = 1 — 1/T. The solution of the above equation is 

m(t) = exp(— et/r) + m cos(u>t — 0), 
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where t n = 2nn/u. The value of Q, at the n-th cycle, can be written as 

Q = Qo exp(--^) = Q exp(-n/r) 
tuj 

Qo is a constant independent of n. The above form shows that Q relaxes exponentially 
with the number of cycles (n) of the oscillating field. The 'relaxation' time T is equal to 
Tuj/2ire. It should be noted here that the dynamic transition occurs at T = 1 in the limit 
h — > . So, for h — > near the dynamic transition point (where the linearisation holds 
good) the behaviour of relaxation time is 

r-e-^cr-T^o^o))- 1 

which shows the power law (exponent is unity) divergence of the 'relaxation' time at the 
dynamic transition point. It may be noted here that the critical slowing down discussed 
above (Q relaxes/ vanishes exponentially) is valid only for Q = 0, not for positive Q. 

(d) Behaviour of dynamic 'specific-heat' near transition point 

The total energy averaged over a cycle of the oscillating magnetic field can be written as 
[7]: E to t = § Hdt/r, where H is the Hamiltonian (equation 2) of the system. The dynamic 
'specific heat' C to t can be defined as the temperature derivative of the total energy E to t- 
Now, if Ctot{= dEtot/dT) is plotted against the temperature, the plot shows (Fig. 4) a very 
sharp peak at the transition point in believe that it will diverge for infinite system. This 
behaviour is similar to that observed in well known equilibrium phase transitions. One 
can detect the dynamic transition and can have an estimate about the transition point 
from the temperature variation of the response like dynamic 'specific heat'. 

In the equilibrium phase transition it is well known that the specific heat is re- 
lated to the fluctuation in energy. What one should expect in the case of this type 
of nonequilibrium transition ? To have a direct answer to this question the fluctuation 
in 'energy' is studied as a function of temperature. The fluctuation in total energy is: 
SE 2 ot = (< E 2 0t > — < E tot > 2 )- Here the symbol < ... > denotes the average over various 
samples (obtained from the different MC samples). If this quantity is plotted against the 
temperature, it would also get sharply peaked at the transition temperature [8]. 

(e) Evidence of diverging length scale 

The evidence of dynamic correlation length of this type of dynamic transition observed 
in the Ising like extended system was also reported [14]. The dynamic susceptibility (in 
2D) is defined as x — L 2 [< Q 2 > — < Q > 2 ], keeping in mind that the fluctuation- 
dissipation theorem holds good [8] for this type of nonequilibrium transition also. This 
is plotted against 1/R = (27r juS)j < r(ho) > where < r(ho) > is average lifetime or 
nucleation time of the system, x was plotted against 1/R for various values of L(— 
64,90,128). The figure shows the peak of x increases in height with increasing system 
size (L). This clearly indicates the finite size effects in x an d implies the existence of a 
divergent length scale associated with the order parameter correlation function near the 
transition point. It is important to note here that this study was done by varying 1/R 
(keeping T and ho fixed) whereas most of the studies on dynamic transition have been 
done by varying the temperature T (keeping uj and ho fixed). However, the results are 
qualitatively invariant under the choice of tunable parameter. This study [14] gives an 
important idea regarding the divergence of 'length scale' at the transition point of this 
dynamic transition. 
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(f) Dynamic phase transition and hysteresis 

(i) Analytic forms of the loop area and the dynamic correlation near the transition point 
In the earlier sections it was discussed that the dynamic transition has a very close resem- 
blance with the magnetic hysteresis. In this section, the relations among the hysteretic 
loss, the dynamic order parameter and the dynamic correlations are dicussed [9]. 
The form of the oscillating magnetic field is 

h{t) = h cos(cjt). (4) 

The dynamic order parameter is defined as 

Q = f m(t)dt, (5) 

which is nothing but the time-averaged magnetisation over a full cycle of the oscillating 
magnetic field. The hysteresis loop area is 

A = — <j> mdh = h cu <j> m{t) sm(ut)dt, (6) 

which corresponds to the energy loss due to the hysteresis. The Dynamic correlation is 
defined as 

C d =< m(t)h(t) > - < m(t) >< h{t) >, 

where < .. > denotes the time average over the full cycle of the oscillating magnetic field. 
Since < h(t) > = 0, one can write 

C d = — i m(t)h(t)dt = — - i m(t) cos(cut)dt. (7) 

Z7T J 2,71 J 

The dynamic correlation has another physical interpretation. For the cooperatively inter- 
acting spin system, this is the negative of the time averaged spin-field interaction energy 
(per spin) in d-dimensions (< Ef >= —T^[d§12i ~i h(t)dt) over a complete cycle of the 
oscillating field. 

In the dynamically disordered (Q = 0) phase and near the transition point, the time 
series of the magnetisation (m(t)) can be approximated as a square wave with a phase lag 
5 with the applied sinusoidal magnetic field. This approximation works well in the low 
temperature region. 

f 1 for < t < r/4 + 5/cu 
m(t) = l -1 for r/4 + 5/cu < t < 3r/4 + 5/uo (8) 
[ 1 for 3t/4 + 5/uj <t <2tt/u, 

where r is the time period of the oscillating field and 5 is the phase lag between mag- 
netisation m(t) and the magnetic field h(t) = h cos(ut). The value of the hysteresis loop 
area can easily be calculated as 

A = 4h sin (5). (9) 

This form of the loop area was also obtained from the assumption that it is approximately 
equal to four times the product of coercive field and remanent magnetization (here the 
remanent magnetisation equal to unity), where the coercive field is identified as h sin(5) 
(the change in field during the phase lag). Considering the same form of the magnetisation 
the dynamic correlation C can also be calculated exactly as 

2h 

C d =^ cos (5). (10) 

7T 
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From the above forms of A and C d it can be written as 



(4/i ) 2 (2/io/vr) 2 




1. 



(11) 



The above relation shows that the loop area A and the dynamic correlation Cd is ellipti- 
cally related to each other. 

The qualitative nature of the dynamic phase boundarycan be realised by considering 
the simplified form of the instantaneous magnetisation in the ordered phase. The dynam- 
ically ordered region (Q ^ 0) can be approximated by considering the following form of 
the magnetization 



In the above simplified approximation, it was considered (since Q ^ 0) that the mag- 
netisation can not jump to the other well, however the value of initial magnetisation is 
reduced by the amount m r . In the real situation it has been observed that this well is 
not fully square (as assumed above in the form of m(t)), it has a cusp like (or parabolic) 
shape. For m r = 2, the above functional form of m{t) will take the form of eqn. 8 and 
one can get the disordered (Q = 0) phase. Taking the above form of magnetisation the 
dynamic order parameter Q can be calculated as Q = (2 — m r )/2. It may be noted that, 
in this simplified approximation the dynamic order parameter Q is independent of phase 
lag S, which is not observed in the real situation (phase lag shows a peak at the transition 
point) [9]. However, this simple picture can anticipate the convex (looking from the ori- 
gin) nature of the dynamic phase boundary. As the temperature increases m r increases 
and it also increases as the field amplitude increases. Since m r increases as ho and T 
increases, in the simplest linearised assumption, one can consider m r is proportional to 
the product of h and T. Demanding, m r = 2 for the dynamic transition (Q — 0), one can 
readily obtain (h )dT d = constant. This equation tells that the dynamic phase boundary 
will be convex having the shape of rectangular hyperbola (xy = c). The convex nature 
of the phase boundary remains invariant even if one assumes that m r is any increasing 
function of both T and h (for example, power law; m r ~ T x Hq, in this particular case 
the equation of the dynamic phase boundary becomes TJ {ho) y d = constant, it is easy to 
see that this gives the convex shape of the dynamic phase boundary). However, this very 
simple assumption can not describe the entire form of the phase boundary accurately, 
particularly near the end points ((h )d = and Td = 0). 

(ii) General relation among Dynamic order parameter, Hysteresis loop area and the Dy- 
namic correlation 

From the usual definitions (given in earlier section) of C d and A, one can write 



m(t) 



1 for < t < r/4 + 5/u 
< l-m r for r/4 + 5/uj < t < 3r/4 + S/u 
1 for 3r/4 + 6/u < t < 2n/u. 



(12) 




where m{oj) = -k= § m(t) exp( 



iut)dt. So, 



C d 



V2^ 



Re {m(u)) 



and 
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A = —h u^/2irlm(m(uj)). 

The general (complex) form of m(u') will be 

m{u') = \m(u')\ exp(i<f>) 

1/2 



A 2 



2tt I h% 



cxp I 



taI1 1 2*C d 



Note that the phase of m{uj') is independent of ho and u. So, Q is related with A and 
Cd as follows 



<9 



m(t)dt 



J du' j> m(u') exp(iu't)dt 



2tvt 



du' 



(^ 2 C 2 A 2 



+ 



h\u' 2 



(13) 



Above equation gives the general relationship among Q, A and Cd- 

It has been observed that the steady response m(t), to a sinusoidally oscillating mag- 
netic field (h(t) = ho cos(c<jt)), is periodic (with phase lag 5) and has the same periodicity 
(r = 2-n/u) of the field. So, one can write m(t) in a Fourier series as 



m(t) = a + a n cos(nut) + b n sm(nut). 

n=l n=l 

From the usual definitions of Q, A and Cd, it is easy to see that 
a = Q, ai = 2C d /h and &! = A/(-Kh ). 
So, one can write 

2Cd A 

m(t) = Q H — - — cos(wt) + H — sin(cjt) + .. 

ft 7r/?,o 



(14) 



(15) 



Keeping only the first harmonic terms (ignoring higher harmonics) one can easily express 
the instantaneous magnetization as 



m(t) = Q + m cos {ut — 5) (16) 

where the amplitude of magnetization is m = [(2Cd/h ) 2 + (A/ (who)) 2 ] 1 / 2 and the phase 
lag is 5 = tan- 1 (^/(2vrC d )). 



(g) Dynamic phase transition and stochastic resonance 

To study the relationship between dynamic transition and the stochastic resonance [19] 
extensive MC simulations were performed [11] in the kinetic Ising model in presence of 
sinusoidally oscillating magnetic field. The frequency is / = 0.001 (kept fixed throughout 
the study). So, one complete cycle of the oscillating field takes r = 1000 MCSS. A time 
series of magnetization m(t) has been generated up to 10 6 MCSS. This time series contains 
10 3 (since r = 1000 MCSS) cycles of the oscillating field. The dynamic order parameter 
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Q has been calculated for each such cycle. So, the statistics (distribution of Q) is based 
on N s = 10 3 different values of Q. 

The statistical distribution P(Q) of dynamic order parameter Q and its temperature 
dependence have been studied [11] closed to the phase boundary to detect the nature[20] 
of the transition. Fig. 5a shows the distributions P(Q) (at fixed value of the field am- 
plitude) for three different values of temperature. Below, the transition the distribution 
shows only two equivalent peaks centered around ±1. Close to the transition point, a 
third peak centered around zero is developed. As the temperature increases slightly, the 
strength of the third peak increases in cost of that of two other (equivalent) peaks. Above 
the transition, only one peak is observed centered around zero. This indicates that the 
transition is first order or discontinuous. 

What is the origin of this kind of first order transition ? To get the answer of this 
question, the time variation of the magnetization m(t) is studied [11] for several cycles 
of the oscillating magnetic field h(t), close to the transition. Sometimes, the system likes 
to stay in the positive well (of the Landau type double well form of the free energy) and 
sometimes it likes to stay in other. It is obvious that the best time for the system to 
switch from one well to the other one, is when the value of the field is optimum ("good 
opportunity") . So, if the system misses one "good opportunity" (first half period of the 
oscillating field) to jump to the other well it has to wait for a new chance (another full 
period of the oscillating field). Consequently, it shows that the residence time (staying 
time in a particular well) can only be nearly equal to an odd integer multiple of the half- 
period (half of the time period of the oscillating field). This leads to two consequences: 

(1) The distribution of the dynamic order parameter Q would be peaked around three 
values (i) Q ~ 0, when the system utilizes "good opportunity" and goes from one well 
to the other (marked 'A' in Fig. 5a), (ii) Q ~ — 1, when the system misses the "good 
opportunity" to go from the negative well to the positive well and it stays for one (or 
more) full period in the negative well (marked 'B' in Fig.5a), (iii) Q ~ +1, when the 
system misses the "good opportunity" to go from the positive well to the negative well 
and spends one (or more) full period in the positive well (marked 'C in Fig. 5a). As a 
result, the distribution of Q would give three distinct peaks centered at +1, -1 and 0. 

(2) The other consequence of this kind of time variation, of magnetization m(t), is 
the "stochastic resonance" [19]. This can be detected from the distribution of residence 
time (the time system spends in a particular well). The distribution (P r ) of residence 
time (r r ) will be peaked multiply around the odd integer multiple of half-period. One 
such distribution is shown in Fig. 6. The distribution shows multiple peaks around the 
odd integer values (500, 1500, 2500, 3500, 4500 and 5500 MCSS) of half-period (r/2=500 
MCSS, of the driving fields). The heights of the peaks decreases exponentially (dotted 
line in Fig. 6 ) with the peak positions. This is the fingerprint of stochastic resonance [19]. 

Figure 5 shows the distributions of the dynamic order parameter Q for three different 
values of the temperature. Here, the field amplitude ho is quite low in comparison with 
that used in the earlier case (Fig. 5a). It shows that, in the ordered region, this gives two 
(equivalent) peaks (Fig. 5b) and as the temperature increases these two peaks come close 
to each other continuously (Fig.5b) and close to the transition (and also above it) (Fig. 5b) 
only one peak (centered around zero) is observed. This feature reveals the continuous or 
second order transition. Hence, it was proved that a tricritical point would exist which 
separates the nature (discontinuous/continuous) of the dynamic transition. However, a 
recent study [15] claims that the existence of TCP [11] is a correct observation but a finite 
size effect. For small system size the distribution of dynamic order parameter Q shows 
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three peaks very close to the transition point. This was observed earlier [11] for small 
system size. However, the distribution of Q has only two peaks near the transition point 
for much larger systems revealing only the continuous nature of the transition. 

(h) Dynamic phase transition for randomly varying field: 

Very recently, an interesting version of this dynamic phase transition has been pre- 
dicted [10] in a ferromagnetic Ising system when the external field on the system varies 
in time stochastically. The long time response (magnetisation) of a kinetic Ising system 
represented by the Hamiltonian (eqn. 2) is studied when the uniform field over the sam- 
ple h(t) varies randomly in time with a uniform distribution bounded between +ho/2 
and —ho/2. In a Monte Carlo simulation study in two dimension, the nature of the re- 
sponse magnetisation (see Fig. 7a and Fig. 7b) is studied with the dynamic order parameter 
Q{= (1/t) f T m(t')dt'; r ^> 1) which is given by the long-time average (over the active 
duration of the magnetic field) of magnetisation. It was found that Q assumes nonzero 
values below a phase boundary line in the ho — T plane, and vanishes continuously at 
the transition boundary (see Fig. 7c). Again, the dynamic symmetry breaking transition 
occurs due to the competing time scales; the relaxation time of the many-body system 
being larger than the switching time of the random field. Such a dynamic transition 
is again a nonequilibrium transition, very similar to that for oscillating fields discussed 
earlier. It may be mentioned that, in a slightly different context, a discrete map version 
of the mean field equation of motion (eqn. 1) with similar stochastically varying field 
h(t) was analysed recently by Hausmann and Rujan [21]. The dynamic transition for a 
randomly varying magnetic field was also studied [10] by solving the mean-field equation 
(eqn. 1) of motion of average magnetisation. 

(i) Athermal dynamic transition in random field Ising model: 

The kind of nonequilibrium dynamic transition discussed so far was assisted by thermal 
fluctuation. An interesting phenomenon, the athermal hysteresis, has been studied [22] 
recently in random field Ising model. Now one may ask, is there any disorder induced dy- 
namic transition observed at T = ? To investigate this the random field Ising model (in 
2D) in presence of oscillating magnetic field was studied [23] at T = by MC simulation. 
A square lattice of linear size L is taken. Each site is labelled by an integer % and carries 
an Ising spin Si (Si = ±1) which interacts with all its nearest neighbours (spins) with 
a ferromagnetic interaction strength J. At each site i, there is a local quenched random 
field hi. The random fields hi are assumed to be independent and identically distributed 
random variables with a rectangular probability distribution P(h,j). The random field hi 
can take any value from —w/2 to +w/2 with the same probability. The width of the 
distribution is w. In addition, there is a uniform (in space) magnetic field h(t) which is 
varying sinusoidally (h(t) = hocos(u>t)) in time. The amplitude and the frequency are 
denoted by h and u respectively. This kind of model is described by the Hamiltonian 

h = - j x; - E hiSi - h(t) £ s h (17) 

<ij> i i 

under the periodic boundary condition. 

The zero-temperature single spin-flip dynamics is specified by the transition rates (W) 

W(Si -> -Si) = T, if AE < and W(Si -> -Si) = 0, otherwise (18) 

where AE is the change in energy due to spin flip. In words the algorithm is: never flip 
the chosen spin if this process would increase the energy and flip otherwise. We have 
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started with all spins are up (Si = +1) as an initial condition and updated the lattice 
sequentially using the above flipping algorithm. One such full scan over the entire lattice 
consists a Monte Carlo step per spin (or MCS). The instantaneous magnetisation (per 
site) m(t) is easily calculated, 

™(0 = ^£^. (19) 

i 

After an initial transient period the intantaneous magnetisation m(t) has been found 
to be stabilised and periodic with the same periodicity of the applied oscillating field. 
For a particular values of h , u and w the dynamic order parameter Q(— \ § m(t)dt) is 
calculated by averaging over 20 different random disorder (quenched) realisations. 

The simulations are performed on a square lattice of linear size L = 100 and a partic- 
ular value of the frequency (uj = 0.01 x 2n) of the oscillating magnetic field. The time is 
measured in units of Monte Carlo steps per spin or MCS and the values of random field 
and the oscillating field are measured in units of interaction strength J. 

It has been observed numerically that, for fixed values of ho and w, the magnetisation 
becomes periodic (in time) with the same periodicity as the applied sinusoidal magnetic 
field. For the smaller values of the quenched disorder (w = 8.0) and the field amplitude 
(ho = 0.5), the magnetisation oscillates asymmetrically about the zero line i.e., the system 
remains in a dynamically symmetry-broken phase. Consequently, the hysteresis (m — h) 
loop resides on the upper half plane formed by h(t) and m(t). So, the time-averaged 
magnetisation over a full cycle of the oscillating field, the dynamic order parameter, is 
nonzero in the symmetry-broken phase. By increasing the field amplitude (for h = 
2.0) keeping w fixed (w = 8.0), it was observed that the system acquires a dynamically 
symmetric phase, i.e., the magnetisation oscillates symmetrically about the zero line. The 
hysteresis loop is also symmetric. Consequently, the value of the dynamic order paramater 
Q is zero in this dynamically disordered (symmetric) phase. 

In the dynamically disordered phase, the dynamic order parameter Q can be kept at 
zero in two ways, either by increasing the random field width w for a fixed field amplitude 
ho or vice versa. So, in the plane formed by the field amplitude (ho) and the width (w) of 
the quenched disorder (random field), one can think of a boundary line, below which Q is 
nonzero and above which it vanishes. Figure 8(a) displays such a phase boundary in the 
ho — w plane obtained by Monte Carlo simulation. The nature (discontinuous/continuous) 
of the transition depends on the value of w and ho on the phase boundary line. The tran- 
sition across the upper part of phase boundary line is discontinuous and it is continuous 
for the rest part of the boundary. A tricritical point on the phase boundary line separates 
these natures. Figure 8(b) demonstrates two typical transitions for two sets of values of 
w and ho lying just in the left and right sides of the tricritical point (TCP). In the case of 
discontinuous transition, the dynamic order parameter Q jumps to a small nonzero value 
and then vanishes continuously. The uncertainty in the location of the TCP on the phase 
boundary are shown by the circle enclosing it. It was not yet checked whether this TCP 
observed here is a finite size effect or not. 



4. Dynamic phase transitions in the Heisenberg ferromagnet: 

(i) Why Heisenberg model ? 

When the nonequilibrium dynamic transition is studied in the Ising model, it has some 
limitations. The Ising model is a special case of general magnetic model [24], for example, 
the Heisenberg model. The Heisenberg model (with ferromagnetic interactions) having 
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uniaxial anisotropy has some general properties which cannot be found in the Ising model. 
But in the limit of infinite anisotropy, the Heisenberg model can be mapped into the Ising 
model. So, the natural expectation is, the Heisenberg model with uniaxial anisotropy 
can be studied to have the detailed and general microscopic view and the results can be 
checked in the limit of infinite anisotropy (which will give the results in the Ising model). 
In this case of dynamic transitions, mainly in the magnetic model system in presence of 
a magnetic field, oscillates sinusoidally in time, the Heisenberg model can serve a better 
role than an Ising model. It would be quite interesting to know the dynamic response of 
uniaxially anisotropic Heisenberg model in presence of a magnetic field, applied in different 
directions. On the other hand, there is another advantage. The results obtained in the 
Ising model are well established [1] . These results can be used to check the results obtained 
in the Heisenberg model by approaching the limit of infinite anisotropy. This prompted 
to study the dynamic transition in the Heisenberg model with uniaxial and single-site 
anisotropy. Recently, the dynamic transition was studied [25] in the uniaxially anisotropic 
ferromagnetic Heisenberg model and it was observed that the dynamic symmetry of the 
order parameter component (along the anisotropy direction) can be broken in presence 
of a magnetic field, applied along the direction which is perpendicular to the direction 
of anisotropy. This transition was named off-axial transition. The transition is found 
to be continuous and the transition temperature increases as the strength of anisotropy 
increases. 

So, the questions naturally arise what would be the difference in the dynamic tran- 
sitions in presence of a field, applied only along the direction of anisotropy ? How the 
symmetry breaking takes place ? What would be the nature (continuous or discontinuous) 
of the transition ? More interestingly, what would happen in the infinitely anisotropic case 
and in the Ising case ? To get the answers of these questions, the researchers studied the 
dynamic phase transition in classical vector spin models. The dynamical phase transition 
in anisotropic XY-ferromagnet in an oscillating magnetic field is studied recently [26] by 
solving the time- dependent Ginzburg-Landau equation. Very recently, it was observed 
theoretically [27] that the symmetry of the vector spin model can be tailored by applying 
oscillating magnetic field. For example [27], depending on the frequency and amplitude 
of the field the Heisenberg ferromagnet can behave like XY-ferromagnet. The dynamic 
transitions in presence of the axial field (i.e., the magnetic field is applied only along 
the direction of anisotropy) and the off-axial field (i.e., the magnetic field is applied only 
along the direction which is perpendicular to the direction of anisotropy) are studied [25] 
by Monte Carlo simulation using Metropolis rate. Also, a comparison between axial and 
off-axial transitions has been made and the results (in the limit of infinite anisotropy) for 
both cases are compared with that observed in the Ising model. By the application of 
polarised magnetic field the multiple dynamic transitions were observed [28] in anisotropic 
ferromagnetic Heisenberg model. 
(ii). The description of the model 

The Hamiltonian of a classical anisotropic (uniaxial and single-site) Heisenberg model 
with nearest neighbour ferromagnetic interaction in the presence of a magnetic field can 
be written as 

H = -JJ2 S t -S J -Dj2(S tz ) 2 -h-J2S t , (20) 

<ij> i i 

where Si[Si X , Si y , Sy represents a classical spin vector of magnitude unity situated at the 
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i-th lattice site. So, Sf x + Sf + Sf z = 1 is an equation of a unit sphere. Classical spin 
means, this spin vector can be oriented in any direction in the vector spin space. J(> 0) 
is the uniform nearest neighbour strength of the ferromagnetic interaction. The factor 
D in the second term is the strength of uniaxial (z here) anisotropy favouring the spin 
to be aligned along the z-axis. The last term is the spin-field interaction term, where 
h[h x , h y , h z ] is the externally applied magnetic field (uniform over the space). When the 
magnetic field is applied only along the a - direction, the magnetic field component h a 
(may be any one of x, y and z) is oscillating sinusoidally in time and can be written as 
h a {t) = /i°cos(cjt), where h° a and uj are the amplitude and angular frequency (cu = 2nf; 
f is frequency) of the oscillating field respectively. Magnetic field \h\ and strength of 
anisotropy D are measured in units of J. The model is defined in a simple cubic lattice 
of linear size L with periodic boundary conditions applied in all the three directions. 

(Hi) The Simulation technique 

The model, described above, has been studied extensively by Monte Carlo simulation using 
the following algorithm [30]. Initial configuration is a random spin configuration. Here, 
the algorithm used, can be described as follows. Two different random numbers r\ and 
T2 (uniformly distributed between -1 and 1) are chosen in such a way that R 2 = (rf + r\ ) 
becomes less than or equal to unity. The set of values of r\ and r 2 , for which R 2 > 1, are 
rejected. Now, u = y/1 — R 2 . Then, S ix = 2wi, S iy = 2ur 2 and S iz — 1 — 2R 2 . 

Starting from an initial random spin configuration (corresponding to high temperature 
configuration) the system is slowly cooled down. At any fixed temperature T (measured 
in the unit of J /Kb) and field amplitude h° a (measured in the unit of J) a lattice site 
% has been chosen randomly (random updating). The value of the spin vector at this 

— * 

randomly chosen site is Si (say). The energy of the system is given by the Hamiltonian 
(equation 1) given above. Now, a test spin vector S[ is chosen randomly (described 
by the algorithm above). For this choice of spin vector at site % the energy will be 
H' = —JJ2<ij> S'i ' Sj — D J2i(S'i Z ) 2 — h-J2i S'i- The change in energy, associated with this 
change in direction of spin vector from Si to S'^ is AH — H' — H. Now, the Monte Carlo 
method [24, 20] will decide how far this change is acceptable. The probability of the change 
is given by Metropolis rate [24, 20] (used here) W(Si -> S-) = Min[l, exp(-AH/K B T)]. 
This probability will be compared with a random number R p (say) between and 1. If R p 
does not exceed W , the move (the change Si — ► 5-) is accepted. In this way the spin vector 
Si is updated. L 3 such random updates of spins, defines one Monte Carlo step per site 
(MCSS) and this is considered as the unit of time in this simulation. The linear frequency 
(/ = u /2n) of the oscillating field is taken 0.001 and was kept constant throughout this 
simulational study. So, 1000 MCSS is required to get one complete cycle of the oscillating 
field and consequently 1000 MCSS becomes the time period (r) of the applied oscillating 
magnetic field. To calculate any macroscopic quantity, like instantaneous magnetisation 
components, the following method was employed. Starting from an initially random 
configuration (which corresponds to a high temperature phase) the system is allowed to 
be stabilised (dynamically) up to 4 x 10 4 MCSS ( i.e., 40 complete cycles of the oscillating 
field) and the averages of various physical quantities are calculated from further 4 x 10 4 
MCSS (i.e., averaged over further 40 cycles of the oscillating field). This is quite important 
to get stable hysteresis loop and it is checked that the number of MCSS mentioned above 
is sufficient to get a stable dynamic phase. Here the total length of this simulation for 
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one fixed temperature T is 8 x 10 4 MCSS (which produces 80 complete cycles of the 
oscillating field). Then the system is slowly cooled down (the value of the temperature T 
has been reduced by small interval) to get the values of the statistical quantities in the 
low temperature ordered phase. Here, the last spin configuration obtained at the previous 
temperature is used as the initial configuration for the new temperature. The linear size 
(L) of the system has been taken equal to 20. The CPU time required for 8 x 10 4 MCSS 
is approximately 25 minutes on an Intel-Pentium-III processor. 

(iv) Off-axial dynamic transition 

The instantaneous magnetisation components (per lattice site) m x = J2iSi/L 3 ,m y = 
J2i Sf I L 3 , m z = J2iSi/L 3 are calculated at each time in presence of magnetic field. The 
time-averaged (over a full cycle of the oscillating field) magnetisation components (the 
dynamic order parameter components) Q x = - § m x dt, Q y = - § m y dt and Q z — - § m z dt 

are calculated by integrating (over the complete cycle of the oscillating field) the 
instantaneous magnetisation components. The total (vector) dynamic order parameter is 

— * 

expressed as Q = xQ x + yQ y + zQ z . 

In this paper, two kinds of dynamic transitions were studied and compared. The 
axial transition means the dynamic order parameter component Q z becomes zero from 
a nonzero value at a finite temperature (the transition temperature) in presence of a 
magnetic field h[0, 0, h z ] applied only along the direction which is parallel to the direction 
of anisotropy. Since the uniaxial anisotropy has been taken along the z-direction, in 
this case, the direction of magnetic field has only nonzero ^-component. The off-axial 
transition [25] is the transition in presence of a magnetic field h[h x , 0, 0] applied only 
along the direction which is perpendicular to the direction of anisotropy. In this case, the 
direction of the magnetic field has only nonzero ^-component. 

In the case of axial transition, the instantaneous magnetisation components are cal- 
culated at any fixed temperature T, strength of anisotropy D and amplitude of axial 
magnetic field h° z . The time-eliminated plot of m z — h z gives the axial hysteresis loop. 
It was observed that at high temperature (T = 2.2) the axial hysteresis loop m z — h z is 
symmetric (symmetric means the loop is distributed about h z axis in such a way that the 
total z-component of magnetization, over a complete cycle of field, vanishes) (fig. 9a). As 
a result Q z = 0. And at low temperature (T = 1.0) the m z — h z loop becomes asymmetric 
(Qz 7^ 0) (fig. 9b). In both cases, the m x — h z and m y — h z loops lie almost along h z axis, 
resulting Q x and Q y equal to zero respectively. Thus a dynamic transition occurs (as the 
temperature decreases) at a certain temperature from a symmetric (Q z = 0; Q = 0) to 
an asymmetric (Q z ^ 0; Q ^ 0) dynamic phase in presence of an axial magnetic field. 

What was observed in the case of off'- axial transition ? Recently studied [25] off- 
axial transition shows similar dynamic transition via breaking the symmetry of m z — h x 
loop in presence of an off-axial field (along perpendicular to the anisotropy direction 
i.e., x-direction) . Here, at high temperature (T = 1.8) the m z — h x loop is symmetric 
(and Q z = 0) and m x — h x loop is also symmetric (Q x = 0) (fig.9c). At some lower 
temperature (T = 0.6), the m z — h x loop becomes asymmetric (Q z = 0) and m x — h x 
loop remains still symmetric (Q x = 0) (fig.9d). In both temperatures Q y = 0. So, here 
also a dynamic transition occurs (as the temperature decreases) at a certain temperature 

— * — * 

from a symmetric (Q z = 0; Q = 0) to an asymmetric (Q z ^ 0; Q ^ 0) dynamic phase 
in presence of an off'- axial magnetic field. Interestingly, it may be noted here that in 
higher temperature the m z — h x loop is 'marginally symmetric' (lies very close to h x 
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axis) rather than a symmetric loop (symmetrically distributed away from and about h x 
line). Strictly speaking, the dynamic transition occurs here from a 'marginally symmetric' 
(loop does not widen up) to an asymmetric phase. One can differenciate the symmetric 
phase from the 'marginally symmetric' phase by considering the loop area of that loop 
whose symmetry breaking is considered in the transition. In the symmetric phase, loop 
is sufficiently widened up resulting in a nonzero loop area. In Fig. 9a, the m z — h z loop 
area is 0.686 (symmetric loop; Q z = 0). But the 'marginally symmetric' loops (m z — 
h x ) have vanishingly small area (0.01) (see Fig.9c) and Q z = 0. It may be noted that, 
in the case of off-axial transition, if the magnetic field is applied along the x-direction 
only (oscillates sinusoidally in time) the m x — h x loop is always symmetric (consequently 
Q x = 0) irrespective of the value of temperature and the strength of anisotropy D (z- 
axis). Similarly, for any field, applied along y-direction only, the m y — h y loop is found 
to be always symmetric (i.e., Q y = 0) irrespective of value of T and D. But in both 
cases, whether the off-axial loops i.e., m z — h x or m z — h y will be symmetric (rather 
'marginally symmetric') or asymmetric that depends upon the values of temperature 
T, anisotropy D and the magnetic field amplitude h° x (or h y ). These results signify 
that without anisotropy the dynamic transition (associated with the dynamic symmetry 
breaking) cannot be observed in the classical Heisenberg model. 

To investigate the dependence of transition temperature on the strength of anisotropy 
(D) in the case of axial transition, the temperature variation of dynamic order parameter 
component Q z was studied for different values of D. Figure 10 shows the temperature 
variation of Q z for different values of D. Here, like the case of off- axial transition [25] the 
transition temperature increases as the strength of anisotropy increases. It is observed 
that the axial transition is discontinuous for lower values of anisotropy strength D (i.e., 
0.5, 2.5 etc.) and it becomes continuous for higher values of D (i.e., 5.0, 15.0 etc.). In the 
Ising limit (D — > oo) the axial transition is also shown in the same figure for D = 400). 
This choice of the value of D{= 400) is not arbitrary. In the case of equilibrium transition 
it was checked by MC simulation that the value of the magnetisation at any temperature 
(in the ferromagnetic region) becomes very close to that obtained in the Ising model, at 
that temperature, if the strength of anisotropy is chosen above 300. 

The temperature variations of dynamic order parameter component Q z in the case 
of off-axial transition was also studied and shown in figure 11 for different values of D. 
This shows that the transition temperature increases as D increases. Here, the transition 
is continuous for all values of strength of anisotropy D. The transition for D = 400 
(D — > oo limit) was compared with that in the case of Ising model. This shows that 
both are continuous and occur at the same point (T 4.5) which is very close to the 
Monte Carlo results of equilibrium ferro-para transition temperature (T c ~ 4.511) [24] in 
3-dimensional Ising model. 

The nonequilibrium dynamical phase transition in the uniaxially anisotropic Heisen- 
berg model, in presence of a magnetic field, which oscillates sinusoidally in time, is studied 
by Monte Carlo simulation. Two cases were studied in this paper, (i) magnetic field, os- 
cillates sinusoidally in time is applied only along the direction of anisotropy, (ii) magnetic 
field applied only along the direction perpendicular to the direction of anisotropy. The 
transition observed in the first case is named axial and that corresponding to the sec- 
ond case is called off-axial. A comparative study between axial and off-axial transition 
is reported in this paper. Three important aspects are considered here, (a) symmetry 
breaking, (b) the order of the transition and (c) the transition in the infinitely anisotropic 
limit. 
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A dynamic symmetry breaking is observed with this dynamic transition. In the case 
of axial transition the dynamic transition occured as the temperature decreases from a 
symmetric to an asymmetric phase, whereas, in off-axial case this symmetry breaking 
takes place from a 'marginally symmetric' to an asymmetric phase. The reason behind it 
is that as follows: in the case of axial transition by the application of axial field (oscillates 
sinusoidally in time) there is a chance that the spin component along the z-direction may 
be reversed in opposite direction which would lead to sufficiently wide and symmetric 
m z — h z loop. But in the case of off-axial transition it is not possible to reverse the 
^-component of spin by applying a field (oscillates sinusoidally in time) perpendicular 
to the direction of uniaxial anisotropy. In this case the value of the ^-component of 
magnetisation m z is almost zero. As a result the m z — h x loop lies on h x = axis and 
hence the loop is marginally symmetric. 

In both cases (axial and off-axial) the transition temperature increases as the strength 
of anisotropy increases provided the amplitude of the applied field remains same. The 
strength of anisotropy tries to align the spin vector along the direction of anisotropy. 
So, as the strength increases it becomes harder to break the symmetry and consequently 
more thermal fluctuation is required to break the symmetry. As a result, the transition 
temperature increases as the strength of anisotropy increases. But the difference is the 
nature of transition. In the case of axial field the transition is discontinuous for lower 
values of anisotropy and it becomes continuous for higher values of anisotropy. The reason 
behind it is, that the axial transition occurs in presence of axial field which reverses 
the z-component of magnetisation. So, in lower values of anisotropy the spin vector 
becomes comparatively more flexible and the transition occurs mechanically in presence 
of axial field at lower temperature and it is discontinuous. As the anisotropy increases 
the effect of axial field (of same value) becomes weak and the transition is driven by 
thermal fluctuations and the transition is continuous. In the case of off-axial transition, 
the off-axial field cannot reverse the z-component of magnetisation. But as the value of 
off-axial field increases, the value of x-component of magnetisation increases at the cost 
of ^-component of magnetisation. The transition is driven by thermal fluctuations and 
continuous. 

What will be the situation in the limit of infinite strength of anisotropy ? In the case of 
axial transition it was observed that the transition temperature for infinitely anisotropic 
Heisenberg model differs from that obtained in an Ising model. Although the equilibrium 
transitions in the infinitely anisotropic Heisenberg model and in the Ising model gives 
the same transition temperature, the nonequilibrium transition temperatures in those 
two cases are not same. Since the magnetic field is applied in z-direction, oscillating 
sinusoidally in time, keeps the system always away from the equilibrium, the system 
does not become an Ising system even in the infinitely anisotropy limit. As a result, the 
dynamic transition temperature in the infinitely anisotropic Heisenberg model cannot be 
same for that obtained in the Ising model. But in the case of off-axial transition, the 
transition temperatures in the infinitely anisotropic Heisenberg model and in the Ising 
model becomes exactly equal. The reason behind it is that as follows: in the case of 
off-axial transition the field is applied perpendicular to the direction of anisotropy. The 
effect of axial field, oscillating sinusoidally in time, has no effect in the infinitely anisotropic 
limit. Though the magnetic field is applied in the x-direction, oscillates sinusoidally in 
time, the infinite anisotropic Heisenberg model becomes an Ising model in statistical and 
thermal equilibrium. Hence, the infinitely anisotropic Heisenberg model in presence of 
off-axial field maps into the Ising model in zero field. That is why the nonequilibrium 
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transition in the infinitely anisotropic Heisenberg model in presence of off-axial field and 
the Ising model (in zero external field) give the same result. 

One important point may be noted here regarding the iclassical dynamics chosen in 
this simulation. Since, the spin component does not commute with the Heisenberg Hamil- 
tonian the spin component has an intrinsic quantum dynamics. Considering this intrinsic 
quantum dynamics there was a study [29] about structure factor and transport properties 
in XY- model. However, in this paper, the motivation is to study the nonequilibrium 
phase transition driven by thermal fluctuations. To study this, one should choose the 
dynamics which arises due to the interaction with thermal bath. Since the objective is 
different, in this paper, the classical dynamics chosen here (which arises solely due to the 
interaction with thermal bath), is Metropolis dynamics. The quantum effect of intrinsic 
spin dynamics is not taken into account. 

(v) Multiple dynamic transition 

In this case the MC study was done [28] for a polarised magnetic field having the 
form: h = ih x + jh y + kh z = ih Qx cos(ut) + kh 0z sm(ut) . One can readily check that 
h x = h 0x cos(ut) and h z = h 0z sm(ujt) yield, after the elimination of time 

h 2 h 2 

w + w = l (21) 

The simulational study is done for a simple cubic lattice of linear size L = 20. The 
total (vector) dynamic order parameter can be expressed as Q = iQ x + jQ y + kQ z . The 
instantaneous energy e(t) = —JJ2<ij> Si • Sj — D J2i(S^) 2 — h -J2iSi is also calculated. 
The time-averaged instantaneous energy is E = § e(t)dt/r. The rate of change of E with 
respect to the temperature T is defined as dynamic specific heat C(= dE/dT) [7]. The 
dynamic specific heat C is calculated from energy E, just by calculating the derivative 
using the three-point central difference formula, given below. 

_ dE _ E(T + ST) - E(T - 5T) 
6 " ~dT ~ 25T (22) 

For the elliptically polarised (equation 21) magnetic field, where the resultant field lies 
in x — z plane, the amplitudes of fields are taken as h 0x = 0.3 and h 0z = 1.0 and the 
frequency / = 0.02. The strength of uniaxial anisotropy is taken D = 0.2. This value of 
D is obtained by rigorously searching to have these interesting results and kept constant 
throughout the study. However, there must be variations in transition points depending 
on the values of D. It is observed that for higher values of D the multiple transition 
phenomenon disappears. The values of field amplitudes and frequency are also obtained 
by searching. 

The temperature variations of the dynamic order parameter components (Q x , Q y , Q z ) 
are studied and the results are depicted in Fig. 12(a). As the system is cooled down, 

— * 

from a high temperature disordered (Q = 0) phase, it was observed that, first the system 
undergoes a transition from dynamically disordered (Q = 0) to a dynamically Y-ordered 
(only Q y ^ 0) phase. This may be called as the first phase (Pi) and the transition tem- 
perature is T cl . This phase can be characterised as Pi. (Q x = 0, Q y ^ 0, Q z = 0). Here, 
the resultant vector of elliptically polarised magnetic field lies in x — z plane and the 
dynamic ordering occurs along y-direction. So, this is clearly an off-axial transition [25]. 
In the case of this type of off-axial transition the dynamical symmetry (in any direction; 
y-direction here) is broken by the application of the magnetic field in the perpendicular 
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direction (lies in the x — z plane here). As the system cools down, it retains this partic- 
ular dynamically ordered phase (Pi) over a considerable range of temperatures. As the 
temperature decreases further, a second transion was observed. Here, the system becomes 
dynamically ordered both in X- and Z-directions at the cost of Y-ordering. In this new 
dynamic phase, P 2 : (Q x 7^ 0, Q y = 0, Q z 7^ 0). In this phase the dynamical ordering is 
planar (lies on x — z plane). The ordering occurs in the same plane on which the field 
vector lies. This transition is axial [25]. This phase may be called the second phase (P 2 ) 
and the transition (from first phase to the second phase) temperature is T c2 . As the tem- 
perature decreases further, the x- and z-ordering increases. At some lower temperature, 
a third transition was observed, from where the x-ordering starts to decrease and only 
^-ordering starts to increase quite rapidly. This third phase can be designated as P 3 : 
(Qx 7^ 0, Q y = 0, Q z 7^ 0). Although the characterisation of P 2 and P 3 , in terms of the 
values of dynamic order parameter components, looks similar there exists an important 
difference between these two phases. In the phase P 2 , both Q x and Q z increases as the 
tempereture decreases but in the phase P 3 , Q x decreases as the temperature decreases 
(see Fig. 12(a)). So these two phases P 2 and P 3 distinctly differ from each other. In this 
phase the dynamical ordering is also axial (along z-axis or anisotropy axis). The system 
continues to increase the dynamical ^-ordering as the temperature decreases further. The 
low temperature phase is only dynamically z-ordered. That means the systems orders 
dynamically (only Q z 7^ 0) along the z-direction (direction of anisotropy) only at very 
low temperatures. Zero temperature dynamic phase (for such polarised field) can be 
characterised as Q x = 0, Q y = and Q z = 1.0. 

To detect the dynamic transitions and to find the transition temperatures the tem- 
perature variation of the energy E is plotted in Fig. 12(b). From this figure it is clear 
that there are three dynamic transitions occur in this case. The transition points are the 
inflexion points in E — T curve. The temperature derivative of the energy E is the dy- 
namic specific heat C. The temperature variation of C is shown in Fig. 12(c). The three 
dynamic transitions are very clearly shown by three peaks of the specific heat plotted 
against the temperature T. From this figure the transition temperatures are calculated 
(from the peak positions of C — T curve). First transition (right peak) occurs around 
T c i = 1.22, the second transition (middle peak) occurs at T c2 = 0.94 and the third (left 
peak) transition occurs around T c3 = 0.86. 

This study was further extended for other values of h 0x keeping other parameters fixed. 
It was found that this three transitions senario disappears for higher values of h 0x . For 
example, for Hq x = 0.9, the second phase P 2 disappears. In this case, the C — T curve 
shows two peaks. It was also observed that for ho x = 0.2, h 0z = 0.2 (keeping all other 
parameter fixed) the system shows single transition and only dynamically orders along 
^-direction. 

In the present study, the external time-dependent magnetic field was taken elliptically 
polarised where the resultant field vector rotates on x — z plane. For the lower values 
of anisotropy and a specific range of the values of field amplitudes the system undergoes 
multiple dynamic phase transitions. Here, three distinct phases are identified. In this 
paper, this observation is just briefly reported. This multiple dynamic phase transition 
in anisotropic Heisenberg ferromagnet in the presence of elliptically polarised field, is 
observed here by Monte Carlo simulation. An alternative method, to check this phe- 
nomenon, may be to use Landau-Lifshitz-Gilbert equation of motion [30] with Langevin 
dynamics. Another important thing should be mentioned here regarding the possible ex- 
planation of multiple dynamic phase transitions (axial and off-axial transitions) observed 
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in the anisotropic Heisenberg model. One possible reason may be the coherent rotation of 
spins. In contrast, the dynamic phase transition in the Ising model can be explained sim- 
ply by spin reversal and nucleation [16]. But to know the responsible mechanism behind 
the multiple dynamic phase transitions, observed in anisotropic Heisenberg ferromagnets 
in the presence of polarised magnetic field, details investigations are required. 

The variation of the dynamic phase boundaries with frequency and the strength of 
anisotropy is quite interesting to be studied. This study also indicates that the system 
will show a very rich phase diagram with multicritical behaviour. The finite size analysis is 
also necessary in order to distinguish the crossover effects from the true phase transitions. 
This requires a huge computational task which will take much time. This work is in 
progress [32] and the details will be reported later. 

In the context of multiple dyanmic transition in anisotropic Heisenberg ferromagnet 
it should be mentioned here that a recent study [31] of anisotropic Heisenberg thin fer- 
romagnetic film shows a double dynamic phase transition for surface and bulk order 
parameter. The Hamiltonian for the classical Heisenberg ferromagnet with a bilinear ex- 
change anisotropy A, in presence of competing surface fields as well as pulsed oscillatory 
fields, was taken [31] as 

H = —JT l< ij > [(l—\)(SiSj+SfSj)+S^Sj)]—T liesU] .{^ e iHi-Si—Tji esuv {^ 

(23) 

where Hi and H D are the static applied surface fields and the time-dependent field H(t) 
was taken to have a pulsed form with 

( _TT 2(k-l)w , < (2fe-l)7T 

Bit) - ~ \ - f <V (24) 

Where ho is the amplitude and u is the angular frequency of the oscillatory external 
field and k is an integer (k = 1,2,3....) representing the number of periods of the pulsed 
oscillatory external field. The model film was taken [31] a simple lattice of size Lx Lx D. 
The system is subject to competing applied surface fields in layers n — 1 and n = D 
of the film with H 1 = hz5n and Hp = —hzSiD. The Monte Carlo study was done [31] 
using Metropolis algorithm for D = 12 and L = 32. They calculated the surface and 
bulk order parameter < Q sur f ace > an d < Q bulk > respectively and studied as a function 
of temperature. From Fig. 13 it is clear that the critical temperature for the dynamic 
phase transition (DPT) in the surface layers is not the same as that for bulk of the film. 
A double DPT was observed [31] for a anisotropic Heisenberg film for competing surface 
fields and pulsed oscillatory fields. 

The off- axial [25] and multiple DPT [28, 31] in anisotropic Heisenberg ferromagnet 
shows quite interesting DPT which was not observed in the Ising model [1] and a rich 
phase diagram is expected here [32]. 

5. Experimental evidences of dynamic phase transitions: 

Several experimental works [33, 34, 35] were performed to investigate the hysteretic 
responses as well as the DPT in ferromagnetic samples. In a recent experiment, Jiang 
et al [35] studied the frequency- dependent hysteresis of epitaxially grown ultrathin (2 to 
6 monolayer thick) Co films on a Cu(OOl) surface at room temperature. The films have 
strong uniaxial magnetisation with two ferromagnetic phases of opposite spin orientations. 
This magnetic anisotropy makes it appropriate to represent the system by an Ising-like 
model. The external magnetic field h(t) on the system was driven sinusoidally in the 
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frequency (/ = u/2n) range 0.1 to 500 Hz and in the amplitude (ho) range 1 to 180 Oe. 
Here of course the time-varying current or the magnetic field induces an eddy current in 
the core, which results in a counter-field reducing the effective magnitude of the applied 
field. The surface magneto-optical Kerr effect technique was used to measure the response 
magnetisation m(t). A typical variation of the loop area A with the driving frequency /, 
at room temperature and at fixed external field amplitude ho, is shown in Fig. 14(a). Fig. 
14(b) shows clearly a signature of DPT and corresponding symmetry breaking. Here, the 
dynamic order parameter Q is plotted against the field amplitude h at a fixed frequency 
and emperature. The inset of Fig. 14(b) shows that for lower values of field amplitude 
i.e., 12.0 Oe (left inset) the phase is dynamically ordered (Q ^ 0) and asymmetric, and 
for higher values of field amplitude i.e., 48.1 Oe, the m — h loop is symmetric and the 
phase is dynamically disordered (Q = 0). This experimental observation supports the 
theoretical results of DPT and dynamical symmetry breaking. However, an experimental 
study of the entire phase boundary is yet to be done. 

6. Summary: 

The DPT in model ferromagnetic systems (Ising and Heisenberg) in the presence 
of sinusoidally oscillating magnetic field is reviewed here. This nonequilibrium DPT is a 
prototype of nonequilibrium phase transition. The kind of nonequilibrium phase transition 
discussed above is observed very recently. Here, only the observations are reviewed. The 
detailed mechanism behind this type of nonequilibrium transition is not yet known clearly. 
In Ising models, the mechanism was tried to be understood in view of nucleation [16]. 
However, in the Heisenberg model the axial, off-axial [25] and very recently observed 
multiple DPT [28] are just observed and the mechanisms responsible for those transition 
are not yet known. The coherent spin rotation [30] may be a possible reason for this. The 
experimental observations of the DPT [35] are only made for Ising like (highly anisotropic) 
systems. Experimental studies are required to observe the special DPT in the Heisenberg 
model. 
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Figure Captions 



Figure 1: Schematic time variation of the response magnetisation m{t) compared to that of the 
oscillating field h(t) for different values of frequency uj and amplitude ho of the oscillating field 
and temperature T of the system. The results are in fact actual Monte Carlo simulation results 
for an Ising model on a square lattice with the values for ho and T as indicated in the Figures. 
The Figures on the right hand side show the corresponding m — h loops. The values for loop 
area A = §mdh and the dynamic order parameter Q are also indicated in these figures. As one 
can see, the first two cases correspond to Q = 0, while the other two correspond to dynamically 
broken symmetric phase (with Q ^ 0). The first figure and the last correspond to A ~ 0, while 
the middle two correspond to nonvanishing A. 

Figure 2: Phase diagrams in the ho-T plane for various values of uj gives the functional form of 
the transition temperature Td{ho,uj) for the dynamic phase transition: Monte Carlo results (a) 
for system sizes L = 100 in d = 2, and (b) for L = 20 in d = 3. Below T^(ho,uj), Q acquires 
a nonzero value in F phase and Q = in P phase. Different symbols denote different phase 
boundary lines corresponding to different frequencies (uj): (□) uj = 0.418, (A) uj = 0.208, (o) uj 
= 0.104 in (a); and (o) uj = 0.418, (□) uj = 0.202, (o) uj = 0.104 in (b). The locations of the 
tricritical points (TCP) are indicated by the circle. The insets show the nature of the transition 
just above (I: ho = 2.2 and 4.4 in (a) and (b) respectively) and below (II: ho = 1.8 and 3.6 in 
(a) and (b) respectively) the tricritical points along the phase boundaries. 

Figure 3: Monte Carlo results of the temperature (T) variation of 'relaxation' time (r) for two 
different values of field amplitudes (ho):the bullet represents ho = 1.0 and the diamond represents 
ho = 0.5. Solid lines show the temperature (T) variations of dynamic order parameter Q. Inset 
shows a typical 'relaxation' of Q plotted against the number of cycles (n). The solid line is the 
best fit exponential form of the data obtained from MC simulation. Here, L = 100, uj = 2-7rx0.04. 

Figure 4: Monte Carlo results of the temperature variations of Ctot for two different values of 
ho- the filled square represents ho = 0.8 and the filled triangle represents ho = 0.4. Solid lines 
represent the temperature variations of Q. Inset shows the temperature variations of E to t for 
two different values of h :(T) h = 0.8 and (II) ho = 0.4. Here, L = 100, uj = 2vr x 0.01. 



Figure 5: (a) The histograms of the normalized distributions of the dynamic order parameter 
Q for different temperatures (T = 0.20J/K B , 0.28J/K B , 0.30J/K B and 0A0J/K B ) and for 
the fixed value of field amplitude ho- All the figures are plotted in the same scales, (b) The 
normalized distributions of the dynamic order parameter Q ( in the 2nd order and close to the 
transition region ) for three different temperatures (T = 1A8J/K B , 1.50 J/ Kb, 1.55J / K b ) and 
fixed field amplitude ho = 0.3J. 



Figure 6: The histogram of normalized (/ P r (T r )dT r = 1) distribution (P r (r r )) of the residence 
time (r r ). The dotted line is the exponential best fit of the envelope of the distribution. 



25 



Figure 7: Dynamic transition due to randomly varing fields in time, (a, b) Typical time 
variation of magnetisation m(t) compared to that of the stochastically varying field h(t) in a 
Monte Carlo study in d = 2, with L = 100, T = 1.7: h = 1.0 for (a) and h = 3.0 for (b). 
(c) The corresponding dynamic transition phase boundary (separating the regions with average 
magnetisations Q zero from nonzero) in the field width (ho) - temperature (T) plane. The data 
points are obtained using both sequential updating (o) and random updating (•) in the Monte 
Carlo simulation. 

Figure 8: (a) The phase boundary (of the dynamic transition) in w — ho plane. The tricritical 
point (TCP) lies within the encircled region. The boundary of the circle is the uncertainty 
associated in locating the TCP, (b) two typical transitions just below and above the tricritical 
point which show the different natures (discontinuous/continuous) of the transitions. 

Figure 9: Symmetry breaking in axial and off-axial transitions. The plot of instantaneous 
magnetization components against the instantaneous field components, (a) m x (t) — h z (t) and 
m z (t) - h z (t) loops for D = 2.5, h° z = 0.5 and T = 2.2, (b) m x (t) - h z (t) and m z (t) - h z (t) 
loops for D = 2.5, h° z = 0.5 and T = 1.0, (c) m x (t) - h x (t) and m z (t) - h x (t) loops for D = 0.5, 
h° x = 0.5 and T = 1.8 and (d) m x (t) - h x (t) and m z (t) - h x (t) loops for D = 0.5, h° = 0.5 and 
T = 0.6. 

Figure 10: The axial dynamic transitions. Temperature (T) variations of dynamic order pa- 
rameter components Q z for different values of anisotropy strength (D) represented by different 
symbols. D = 0.5(O), D = 2.5(+), D = 5.0(a), D = 15.0(x) and D = 400.0(A). In all these 
cases for the axial transitions h Q z = 0.5. The data for the temperature variation of dynamic order 
parameter in the Ising model (for h\ = 0.5 and / = 0.001) are represented by *. Continuous 
lines in all cases are just connecting the data points. 

Figure 11: The off-axial dynamic transitions. Temperature (T) variations of dynamic order 
parameter components Q z for different values of anisotropy strength (D) represented by different 
symbols. D = 0.5(O), D = 2.5(+), D = 5.0(a), D = 15.0(x) and D = 400.0(A). In all these 
cases for off — axial transitions h° x = 0.5. The data for the zero-field ferro-para equilibrium 
Ising transition are represented by *. Continuous lines in all cases are just connecting the data 
points. 

Figure 12: (a) The temperature variations of the components of dynamic order parameters. 
Different symbols represent different components. Q x (diamond), Q y (circle) and Q z (bullet). 
This diagram is for D = 0.2 and for elliptically polarised field where ho x = 0.3 and ho z = 1.0. 
The size of the errorbars of Q x , Q y and Q z close to the transition points is of the order of 
0.02 and that at low temperature (e.g., below T = 0.5) is around 0.003. (b) The temperature 
variation of the dynamic energy (E) for D = 0.2, ho x = 0.3 and ho z = 1.0. The vertical arrows 
represent the transition points, (c) The temperature variation of dynamic specific heat (C = ^) 
for D = 0.2, ho x = 0.3 and ho z = 1.0. Vertical arrows show the peaks and the transition points. 

Figure 13: Surface order parameter < Q sur f ace > an d bulk order parameter < Q bulk > for the 
film, plotted as a function of temperature (T*) for the value of pulsed oscillatory field Hq = 0.3. 
[After Jang et al, Phys. Rev. B 67 (2003) 094411] 
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Figure 14: Experimental results for the dynamic hysteresis loop area A and the dynamic order 
parameter Q [35] . (a) The results for the loop area A as a function of frequency / is plotted at 
a fixed ac current of 0.4 Amp. The direction of the magnetic field is parallel to the film plane. 
The insets show plots of m — h loops for the following particular values of the field amplitudes 
h$: (i) /io = 48.0 Oe (top inset) and (ii) ho = 63.0 Oe (bottom inset), (b) The dynamic order 
parameter Q, i.e, the average magnetisation over a cycle, is plotted against the field amplitude 
at a fixed frequency / = 4 Hz. The insets show plots of m — h loops for the following particular 
values of the field amplitudes ho- (i) ho = 48.1 Oe (right inset) and (ii) ho = 12.0 Oe (left inset). 
[After Q. Jiang et al. Phys. Rev. B. 52 (1995) 14911] 
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